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Abstract. We report the discovery of a new planet candidate orbiting the subgiant star HD 118203 with a period of 
P = 6.1335 days. The best Keplerian solution yields an eccentricity e= 0.31 and a minimum mass m2 sin / = 2. 1 Mjup for the 
planet. This star has been observed with the ELODIE fiber-fed spectrograph as one of the targets in our planet-search pro- 
gramme biased toward high-metallicity stars, on-going since March 2004 at the Haute-Provence Observatory. An analysis of 
the spectroscopic line profiles using line bisectors revealed no correlation between the radial velocities and the line-bisector 
orientations, indicating that the periodic radial-velocity signal is best explained by the presence of a planet-mass companion. 
A linear trend is observed in the residuals around the orbital solution that could be explained by the presence of a second 
companion in a longer-period orbit. We also present here our orbital solution for another slightly evolved star in our metal-rich 
sample, HD 149143, recently proposed to host a 4-d period Hot Jupiter by the N2K consortium. Our solution yields a period 
^=4.09 days, a marginally significant eccentricity e = 0.08 and a planetary minimum mass of 1.36 Mj„p. We checked that the 
shape of the spectral lines does not vary for this star as well. 
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1. Introduction 

Stars hosting planets are significantly metal rich in compari- 
son to field stars in the solar neighbourhood (Gonzalez 1997; 
Santos etal. 2001, 2003, 2005; Fischer & Valenti 2005). These 
authors have also shown that the probability of hosting a giant 
planet is a strongly rising function of the star metal content. 
According to their estimate, we can expect that up to 25 - 30% 
of the more metal-rich stars ([Fe/H] >0.2-0.3) host a giant 
planet. 

On the basis of this argument, a new survey was started in 
March 2004 at the Haute-Provence Observatory with the high- 
precision ELODIE fiber-fed echelle spectrograph. The main 
idea of this new programme is to bias our target sample to- 
wards high-metallicity stars, much more likely to host planets. 
This will strongly increase our probability of finding new plan- 
ets in a sample of yet non-observed stars. The survey uses the 
cross-coiTelation technique for the radial- velocity and metallic- 
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telescope, both at the Observatoire de Haute Provence, France 



ity estimates. The programme mainly targets giant planets with 
short periods (Hot Jupiters). They are the ideal candidates in the 
search for photometric transits. From this survey, we present 
here a new short-period planet candidate (P = 6.1335 d) with 
an eccentric orbit around the star HD 11 8203. 

A similar planet-search programme was simultaneously 
started by the N2K consortium (Fischer et al. 2004) aim- 
ing at the detection of short-period planets orbiting metal-rich 
stars. Fischer (2005, also in Fischer et al. 2005) recently an- 
nounced two new Hot- Jupiter detections around HD 149143 
and HD 109749. HD 149143 is amongst the stars already ob- 
served in our sample and we present here our orbital descrip- 
tion of the system. 

Together, the N2K and ELODIE metaUicity-biased planet- 
search programmes have detected five new Hot Jupiters in less 
than one year. One of them, HD 149026, is transiting in front of 
its parent star (Sato et al. 2005) allowing for the determination 
of the planet radius and mean density. The planet is found to 
have an unexpectedly large core. This result clearly illustrates 
the importance of such programmes for our understanding of 
planet interiors. However, when examining possible statistical 
trends between orbital and stellar parameters to derive con- 
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straints for planet formation models, we have to keep in mind 
the built-in bias of this subsample of exoplanets. In particular 
these planets must be removed when considering correlations 
with the star metallicity. 

The sample selection and observations of the new ELODIE 
programme are described in Sect. 2. Stellar parameters, radial- 
velocity measurements of HD 1 18203 as well as the orbital so- 
lution derived for the new Hot Jupiter candidate are presented 
in Sect. 3. This section also provides information on stellar ac- 
tivity and in particular the results of the bisector analysis for 
the star. Section 4 reports our results for HD 149 143 and con- 
clusions are presented in Sect. 5. 



2. Sample and observational strategy 

The star sample was selected from the HIPPARCOS catalogue 

(ESA 1997). Initially, the bright northern dwarf stars (V < 8.5) 
of colours compatible with spectral types ranging from F8 to 
MO (0.45 < B - y < 1.4) were considered. Then: 

i) Evolved stars were removed from the list, using a criteria 
of proximity to the main sequence in the HR diagram (2 mag). 

ii) All stars already members of the known main pro- 
grammes of search for extra-solar planets in the northern hemi- 
sphere were removed. 

iii) We eliminated stars known as binaries from previous 
CORAVEL measurements. 

iv) We also eliminated stars with high rotational velocities 
in order to avoid active stars in the sample. In practice, stars 
with V sin / > 5 kms ' were excluded from the sample. 

The final selection comprises 1061 stars. They lie at distances 
up to lOOpc from the Sun. 

Observations were conducted with the high-precision 
ELODIE fiber-fed echelle spectrograph (Baranne et al. 1996) 
mounted on the 1.93-m telescope at the Haute-Provence 
Observatory (France). The spectra have a resolution (A/AA) 
of about 42,000. Typical signal-to-noise ratios obtained in 
20 minute exposures range from ~20 to 100 for the pro- 
gramme stars, corresponding to photon-noise errors between 
5 and 20 ms"' on individual measurements. The data reduction 
is performed on-line during the observations by the automatic 
reduction software (see Baranne et al. 1996, for a detailed 
description). 

After a single radial- velocity measurement, it is possible to 
obtain a very good estimate of the star metallicity by measur- 
ing the surface of the cross correlation function (CCF) of the 
ELODIE spectra (Santos et al. 2002; Naef 2003). The typi- 
cal uncertainty of the resulting metallicities is ~0.05 dex com- 
pared to values obtained from a high-resolution spectroscopic 
analysis. The expected percentage of > 10-30% of giant plan- 
ets orbiting stars with [Fe/H] > 0.10 (Santos et al. 2001, 2004; 
Fischer & Valenti 2005) should lead to the discovery of a few 
tens of new planets in our programme, about ten of them in 
short-period orbits. 



HD 118203 ELODIE 




Fig. 1. /i3968.5 A Ca II H absorption line region of the summed 
ELODIE spectra for HD 118203. No clear emission feature is ob- 
served. For clarity, spectral features due to pollution by the thorium 
spectrum have been removed (e.g. around /13966.5 A). 



3. A short-period planet orbiting HD 118203 

3.1. Stellar characteristics of HD 1 18203 

HD 118203 (HIP66192) is Usted in the HIPPARCOS catalogue 
(ESA 1997) as a KO dwarf in the northern hemisphere, with a 
visual magnitude V = 8.05, a colour index B-V = 0.699, and an 
astrometric parallax ;r = 1 1 .29 ± 0.82 mas, setting a distance of 
88.6 pc from the Sun. The corresponding absolute magnitude 
My = 3.31 is too high for a KO dwarf, indicating that the star 
is already sUghtly evolved and is in a subgiant stage. This is 
confirmed by the astrometric surface gravity log ^ = 3.87 dex 
estimated from the eff'ective temperature and stellar mass (see 
below) and Hipparcos precise parallax (Santos et al. 2004). 

From calibrations of the width and the surface of the 
ELODIE cross-correlation functions (Santos et al. 2002; 
Naef 2003), we have estimated a projected rotation velocity 
vsin; =4.7kms~^ and a metallicity [Fe/H] = 0.10 for the star. 
Comparing these observable stellar characteristics with the 
stellar evolution models of Girardi et al. (2002), we can infer 
the following intrinsic properties for HD 11 8203: a coarse esti- 
mate of the effective temperature Teff = 5600+ 150K, a mass 
M* = 1.23 + 0.03 Mo and an age of 4.6 ± 0.8 Gyr. A similar 
value of the effective temperature T^g = 5695 + 50 K is also de- 
rived from the calibration in Santos et al. (2004) using the 
color index and metalUcity of the star. The stellar parameters 
are gathered in Table 1 . 

A check for stellar chromospheric activity was also per- 
formed by looking at the Call H absorption line in the spectra. 
Figure 1 shows the co-added ELODIE spectra where the clear 
absence of an emission feature indicates a low activity level, 
as expected for slightly evolved subgiants (Wright 2004). The 
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Table 1. Observed and estimated parameters for HD 1 18203. (See text 
for references on the quoted values.) 



Spectral Type 


KO 




y 


8.05 




B-V 


0.699 




n 


11.29 ±0.82 


[mas] 


My 


3.31 






5600 ± 150 


[K] 




1.23 ± 0.03 


Mo 


age 


4.6 + 0.8 


Gyr 


logg 


-3.87 


cgs 


[Fe/H] 


0.10 ± 0.05 




vsin (' 


4.7 


[kms-'] 



activity-inducedradial-velocityjitteristhus also expected to be 
low for this moderately rotating star. 

3.2. Orbital solution for HD 1 18203 b 

ELODIE observations of HD 118203 were conducted from 
May 2004 (JD = 2453 151) to July 2005 (JD = 2453 553), 
stimulated by the difference in radial velocity found among the 
first three observations. A set of 43 precise radial- velocity mea- 
surements were then gathered. They are provided in Table 2. 

A clear 6.1-d periodic variation is seen in this data set. 
The residuals around a Keplerian fit with this period are 
however unexpectedly large (>40ms ') and shows a clear 
additional radial-velocity drift as a function of the Julian 
date. A simultaneous Keplerian -i- linear-drift adjustment 
then provides a very satisfactory solution for the system 
with a period f = 6.1335 + 0.0006 days and an eccentricity 
e = 0.309 ±0.014. The slope of the linear drift is found to be 
49.7 + 5.7 ms^'yr"'. It is most probably accounted for by the 
presence of a second companion in the system, on a longer- 
period orbit. This could explain the significantly non-zero value 
found for the eccentricity at such a short period'. 

A plot with the last-season radial- velocity measurements of 
HD 118203 is shown in Fig. 2, together with the derived solu- 
tion. The residuals around the solution are displayed as well in 
the bottom panel of the figure. The phase-folded radial-velocity 
curve with the complete set of data points corrected for the 
observed linear drift is displayed in Fig. 3. The weighted rms 
around the solution is cr(O-C) = 18.1 ms slightly larger than 
the individual photon-noise errors (~15 ms~'). The orbital ele- 
ments are listed in Table 3 with the inferred planetary parame- 
ters. 

3.3. Low chromospheric activity for HD 1 18203 

Physical events in the stellar atmosphere, like the presence of 
spots on the stellar surface, can change the observed spectral- 
Une profiles and induce a transient periodic radial- velocity sig- 
nal similar to the one expected from the presence of a planet. 

' The period is however close to (higher than?) the observed cir- 
cularisation period determined from the sample of known exoplanets 
around 6 days (Halbwachs et al. 2004). 



Table 2. ELODIE radial velocities of HD 118203. All data are rela- 
tive to the solar system barycentre. Given uncertainties correspond to 
photon-noise errors. 



JD-2400000 


RV 


Uncertainty 


[days] 


[kms '] 


[kms '] 


53151.4302 


-29.263 


0.013 


53153.4235 


-29.231 


0.014 


53155.3946 


-29.640 


0.011 


53217.3409 


-29.523 


0.011 


53218.3429 


-29.320 


0.011 


53219.3450 


-29.243 


0.013 


53220.3498 


-29.215 


0.012 


53224.3411 


-29.300 


0.009 


53226.3539 


-29.152 


0.017 


53227.3547 


-29.232 


0.026 


53228.3326 


-29.533 


0.012 


53229.3231 


-29.549 


O.OIl 


53230.3545 


-29.338 


0.025 


53231.3597 


-29.257 


0.013 


53232.3270 


-29.209 


0.018 


53233.3267 


-29.261 


0.009 


53234.3264 


-29.474 


0.010 


53392.6777 


-29.201 


0.014 


53394.6422 


-29.598 


0.027 


53396.6787 


-29.257 


0.015 


53397.6810 


-29.153 


0.015 


53398.6756 


-29.202 


0.015 


53399.6799 


-29.369 


0.014 


53421.6764 


-29.187 


0.024 


53422.6439 


-29.189 


0.017 


53424.6543 


-29.556 


0.013 


53425.6554 


-29.543 


0.014 


53426.6586 


-29.351 


0.012 


53428.6526 


-29.190 


0.025 


53429.5865 


-29.166 


0.028 


53430.5650 


-29.446 


0.016 


53431.6340 


-29.570 


0.012 


53486.5332 


-29.565 


0.012 


53487.4809 


-29.386 


0.011 


53488.4736 


-29.260 


0.012 


53489.5168 


-29.184 


0.017 


53490.4763 


-29.185 


0.010 


53491.4789 


-29.298 


0.011 


53492.4688 


-29.594 


0.011 


53517.4910 


-29.547 


0.011 


53519.4450 


-29.212 


0.011 


53520.3896 


-29.152 


0.011 


53553.3845 


-29.463 


0.025 



The bisector analysis is one of the best tools to discriminate be- 
tween radial- velocity variations due to changes in the spectral- 
Une shapes and variations due to the real orbital motion of the 
star (see Queloz et al. 2001, for a description of this method). 

The bisector inverse slope (BIS value) computed from the 
HD 118203 spectra are plotted in Fig. 4 in comparison to the 
corresponding radial- velocities. If the hne-shape and the radial- 
velocity variations share the same origin, the BIS and radial 
velocities are expected to be correlated. This is not the case. 
Furthermore, no coherent signal is observed when the BIS are 
phased with the orbital period of 6.1335 days. 
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Fig. 2. Top: last 5-months radial-velocity measurements of 
HD 118203 superimposed on the best Keplerian + linear drift 
solution indicative of a second longer-period companion in the 
system. Bottom: residuals around the solution. Error bars represent 
the photon-noise errors. 

Table 3. ELODIE best Keplerian orbital solution obtained for 
HD 118203 as well as the inferred planetary parameters. For eccen- 
tric orbits T is defined as the time of the peri-astron passage. 



p 


6.1335 +0.0006 


[days] 


T 


53394.23 ± 0.03 


[JD-2400000] 


e 


0.309 ± 0.014 




V 


-29.387 ± 0.006 


[km s-'] 


CO 


155.7 ±2.4 


[deg] 


K 


217 + 3 


[m s-'] 


N 


43 




o-(O-C) 


18.1 


[m s-'] 


Linear drift 


49.7 + 5.7 


[m s^'yr"'] 


flisini 


11.61 


[10-5 AU] 


f{m) 


5.55 


[lO-'Mo] 


m2sin i 


2.13 


[Mj„p] 


a 


0.07 


[AU] 



IVIoreover, data gathered during the still negative search 
for a photometric transit (presented elsewhere) show that the 
star is constant at a 0.0047 mag level and thus further support 
its low-activity level. Consequently, the explanation of the ob- 
served radial-velocity variations due to stellar activity can be 
discarded. 



4. A Hot Jupiter around HD 149143 

During the preparation of this publication, another new Hot 
Jupiter candidate was identified in our programme, orbiting the 
star HD 149143. While gathering more points to characterize 
the system, we learned that the planet had just been announced 



Fig. 3. Top: phase-folded ELODIE radial- velocity measurements for 
the star HD 1 18203 after correction for the linear radial- velocity drift 
observed in the data. Bottom: corresponding residuals around the so- 
lution. Error bars represent the photon-noise errors. 
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Fig. 4. Radial velocities (RV, upper panel) and inverse bisector slope 
(BIS, middle panel) phased with the orbital period P = 6.1335 d for 
HD 118203. The independence of the two quantities is shown by the 
radial velocity vs. BIS plot (bottom panel). 
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Fig. 6. Phase-folded ELODIE radial velocities of HD 149143 super- 
imposed on the derived quasi-circular Keplerian solution. Residuals 
around the solution are displayed in the bottom panel. Error bars rep- 
resent the photon-noise errors. 



Fig. 5. ELODIE radial velocities of HD 149143 superimposed on the 
derived Keplerian quasi-circular model. Residuals around the solution 
are displayed in the bottom panel. Error bars represent the photon- 
noise errors. 

Table 4. Observed and estimated parameters for HD 149143. (See text 
for references on the quoted values.) 



Spectral Type 


GO 




V 


7.9 




B-V 


0.68 




n 


15.75 + 1.07 


[mas] 


Mv 


3.88 






5730 ± 150 


[K] 




1.1 ±0.1 


Mo 


age 


7.6 ± 1.2 


Gyr 


logg 


4.10 


cgs 


[Fe/H] 


0.20 ± 0.05 




vsin ; 


3.9 


[kms-'l 



by the N2K consortium (Fischer 2005; Fischer et al. 2005). 
We thus present here our ELODIE solution for the system, rely- 
ing on the N2K photometric measurements for the transit non- 
detection. 

4.1. HD 149143: stellar characteristics 

HD 149143 is Hsted in the Hipparcos catalogue (HIP81022) 
with a GO spectral type, a visual magnitude V = 7.90, and a 
colour index B- V = 0.68. From calibrations of the width and 
the surface of the ELODIE cross-correlation functions we have 
estimated a projected rotation velocity vsin / = 3.9 kms"' and a 
metallicity [Fe/H] =0.2. 

The measured Hipparcos parallax of HD 149143 is 
7r= 15.75 + 1.07 mas setting the star 63.5 pc away from the 
Sun. The inferred absolute magnitude is then My = 3.88, also 
slightly overluminous for a metallic GO dwarf. This is con- 



firmed by the astrometric estimate of log g-4.l0 dex. From the 
models of Girardi et al. (2002) we obtain an effective temper- 
ature Teff = 5730 +150 K, a primary mass Mi = 1.1 +0.1 Mq 
and an age of 7.6 + 1.2 Gyr. A calibration of the stellar effec- 
tive temperature based on the stellar colour and metallicity 
(Santos et al. 2004, only slightly sensitive to the gravity of 
the star) leads then to Teff = 5790 + 50 K, compatible with the 
value given above. The stellar parameters are given in Table 4. 

4.2. Orbital solution for HD 149143 b 

Although the number of ELODIE observations is small for 
this star (Table 5), thanks to the large amplitude of the radial- 
velocity variation (/T = 163 ms"'), we can easily derive a quasi- 
circular Keplerian solution (e - 0.08 + 0.04) with a period of 
4.088 days, in agreement with the Fischer (2005) and Fischer 
et al. (2005) announcement. The eccentricity of the orbit is 
small and marginally significant, we thus also provide the cor- 
responding circular solution fixing e=0 (Table 6). The derived 
orbital elements coupled with the above estimate of the primary 
mass of 1.1 Mo lead to a minimum mass m2 sin/= 1.36 Mjup 
and a separation of 0.052 AU for the planetary companion. 

Figures 5 and 6 present the temporal and phased ELODIE 
radial-velocity measurements superimposed on the derived 
Keplerian model as well as the residuals around the solu- 
tion. The inferred orbital and planetary parameters are given 
in Table 6. 

As for the case of HD 1 18203, the bisector inverse slope of 
the cross-coiTelation function has been calculated. No coiTela- 
tion is found between these bisector slopes and the radial ve- 
locities, excluding activity-induced variations of the shape of 
the spectral lines as the source of the radial-velocity variations. 

5. Summary and concluding remarks 

We have presented the characteristics of a new planet candidate 
in orbit around the subgiant star HD 118203, detected by the 
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Table 5. ELODIE radial velocities of HD 149143. All data are relative 
to the solar system barycentre. 



JD-2400000 


RV 


Uncertainty 


[days] 


[kms '] 


[kms '] 


53550.4467 


12.226 


0.013 


53554.4224 


12.174 


0.010 


53587.3709 


12.232 


0.014 


53589.3393 


11.880 


0.013 


53592.3736 


12.102 


0.012 


53611.3172 


12.121 


0.010 


53613.3130 


11.992 


0.010 


53614.3141 


11.928 


0.010 



Table 6. ELODIE derived Keplerian orbital solution obtained for 
HD 149143 as well as the inferred planetary parameters. Both the 
quasi-circular and circular solutions are given because the derived ec- 
centricity is only marginally significant. For the eccentric orbit T is 
defined as the time of the peri-astron passage whereas for the circular 
orbit T indicates the maximum of radial velocities. 



e free e = 

f 4.088 ± 0.006 4.089 ± 0.006 [days] 

T 53588.0 + 0.4 53587.4 ±1.0 [JD-2400000] 

e 0.08 ± 0.04 0. 

V 12.059 ±0.004 12.056 ±0.003 [kms '] 

Oi 42 ± 35 0. [deg] 

K 163 ±8 156 ±6 [ms"'] 

o-(O-C) 13.3 16.0 [ms-'] 



flismi 6.105 5.847 [10"^ AU] 

/(w) 1.82 1.59 [10"Mo] 

mzsini 1.36 1.30 [Mj^p] 

a 0.052 0.052 [AU] 



new ELODIE planet-search programme biased towards metal- 
rich stars. The planet is in a rather eccentric orbit (e = 0.31), 
with a period of P = 6.1335 days, and is close to its parent star 
(a = 0.06AU). 

An additional trend of the radial-velocity measurements 
increasing as a function of JuUan date with a slope of 
49.7 ms-Vr~^ is observed, suggesting the presence of a sec- 
ond companion around the main star on a longer-period orbit. 
This additional companion could explain the somewhat high 
eccentricity of the orbit of this new planet at short period. 

We have also reported the ELODIE solution for another Hot 
Jupiter in our programme, HD 149143, recently announced by 
the similar-goal N2K project. The derived best Keplerian quasi- 
circular solution presents a period of 4.09 days and the inferred 
planetary mass is 1.36Mjup. 

By selection these planets are orbiting metal-rich stars. 
They increase to five the number of Hot Jupiters detected in 
less than one year by dedicated metallicity-biased programmes. 
One of these candidates, HD 149026 (Sato et al. 2005), transits 
in front of its parent star and thus allows the determination of 
its radius and mean density when combining photometric and 
radial- velocity measurements. This demonstrates the efficiency 
of such approaches to find candidates suitable for constraining 
planet-interior models. However, the built-in biases of the sam- 



ple have to be kept in mind when examining possible statistical 
relations between the star metalUcity and other orbital or stellar 
parameters. 
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